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Abstract 

In e + e“ collisions recorded using the CLEO II.V detector we have studied the Cabibbo sup¬ 
pressed decay of D° —> 7r + 7r _ 7r° with the initial flavor of the D° tagged by the decay D* + —> D°tt + . 
We use the Dalitz-plot analysis technique to measure the resonant substructure in this final state 
and observe pir and non-resonant contributions by fitting for their amplitudes and relative phases. 
We describe the 7T7 t S'-wave with a iC-matrix formalism and limit this contribution to the rate to 
be < 2.5% @ 95% confidence level, in contrast to the large rate observed in D + —> 7r + 7r - 7r + decay. 
Using the amplitudes and phases from this analysis, we calculate an integrated CP asymmetry of 
O.Olloio? ± 0-05. 


2 


The light scalar meson sector is an enduring puzzle in QCD JJ. Specifically the isospin 
zero, J PC = 0 ++ mesons are complex both theoretically and experimentally. The singly- 
Cabibbo suppressed decays of D mesons are an excellent laboratory to test this sector. The 
final state consists of only u and d quarks and antiquarks so there is sufficient energy to 
cover most of the range of interest to light quark binding, and the initial state is simple 
with little impact on the final state. A better understanding of final state interactions in 
exclusive weak decays is needed in order to model rates, explain interesting phenomena such 
as mixing '2| and elucidate the origin of CP violation in the B sector 

Weak decays of D mesons are expected to be dominated by resonant two body decays jij, 
Si HI- The well established Dalitz-plot analysis technique can be used to explore the 
resonant substructure which should be rich in isospin zero mesons. Recently the FOCUS 
collaboration studied the Dalitz plot D + —> 7r + 7r - 7r H 
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and observed a large ini S -wave 
contribution using a Jl-matrix formulation to describe to 0 ++ resonance structures. This 
letter describes a similar Dalitz-plot analysis of D° —> 7r + 7r _ 7r° at CLEO in which nn S- 


wave contributions are also expected. We have searched for such contributions using the K- 
matrix formulation following what has been done by FOCUS and alternatively for resonance 
contributions from the scalar cr(500) and /o(980) mesons. We see no evidence for any mr 
S -wave contribution in D° —> 7r + 7r _ 7r°, and fully describe the Dalitz plot with contributions 
from p resonances. 

Standard Model (SM) predictions for the rate of CP violation in almost all charm meson 
decay modes is 0( 10 -6 ). However, for some singly Cabibbo suppressed decays of D mesons 
such as D° —> 7T + 7r - 7r°, the SM predictions for the rate of CP violation are as large as 


0.1% 8. |!)|. due to interference between tree and penguin processes. 

Previous investigations El of this decay were limited by statistics and did not search for 
CP violation nor study the resonant substructure. 

This analysis uses an integrated luminosity of 9.0 fb -1 of e + e _ collisions at a/s ~ 10 GeV 
provided by the Cornell Electron Storage Ring (CESR). The data were taken with the 
CLEO II.V detector QQ. 

We reconstruct candidates for the decay sequence D* + — >7r p D°, P°-> 7r + 7r _ 7r°. Charge 
conjugation is implied throughout this Letter. The charge of the slow pion (jr* or 7r~) 
identifies the charm state at t = 0 as either D° or D°. To reduce background, we require the 


D* + momentum to exceed 70% of its maximum value yjE^ eam — M|,* + . The 7r° candidates 
are reconstructed from all pairs of electromagnetic showers that are not associated with 
charged tracks. To reduce the number of fake 7T° from random shower combinations we 
require that each shower has energy greater than 100 MeV and be in the barrel region of our 
detector. The two photon invariant mass is required to be 120 < M 77 < 150 MeV/c 2 . To 
improve the mass resolution, the invariant mass is constrained to the known n° mass and 
we require the y 2 of this fit to be < 100. We exploit the precision tracking of the silicon 
vertex detector [l3[ by refitting the 7r ± tracks with a requirement that they form a common 
vertex in three dimensions. We use the trajectory of the 7r + 7r _ 7r° system and the position of 
the CESR luminous region to obtain the D° production point. We then refit the 7T+ track 
with a requirement that the trajectory intersect the D° production point. 

We reconstruct the energy released in the D* + — >n p D° decay as Q = M*—M—m n , where 
M* is the reconstructed mass of the 7r/7r + 7r _ 7r° system, M is the reconstructed mass of the 
7T + 7r - 7r° system, and m r is the charged pion mass. The addition of the D° production point 
to the 7r/ trajectory improves the resolution on Q by a factor of two. The distributions of 
Q and M for our data are shown in Fig. |T| We fit the M and Q distributions separately 
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to a double Gaussian plus a background shape and find an average background fraction of 
(18.6 ±3.6)%. We select 1917 candidates within 650 keV of the nominal value of Q , denoted 
as Q, and within 44 MeV/c 2 of the nominal value of M, both as measured in this analysis. 

The efficiency for the selection described above is not uniform across the Dalitz-plot dis¬ 
tribution (m 2 +7r _, m 2 +?r0 ). We study the efficiency with a GEANT |l4j based simulation of 
the detector with a luminosity corresponding to more than twenty times our data sample. 
To measure the variation in efficiency over the Dalitz plot, we generate signal Monte Carlo 
uniformly populating the allowed phase space. We observe deviations from the uniform dis¬ 
tribution due to momentum dependent n° reconstruction efficiency and inefficiencies near 
the edge of phase space. The average reconstruction efficiency is ~3.7% but increases (de¬ 
creases) to ~4.5% (~0.5%) for decays with high (low) momentum 7r° mesons. We fit the 
efficiency to a two dimensional cubic polynomial in (m 2 +7r _, m 2 +7r0 )- 

Figure |H shows that the background is significant. To construct a model of the background 
shape, we consider events in the data in sidebands 3 < Q — Q < 6 MeV and 3 < Q — 
Q < 15 MeV within the M signal region defined above. There are 2711 events in this 
selection, about eight times the amount of background we estimate from the signal region. 
The backgroimd is dominated by random combinations of unrelated tracks and showers. 
Although the background includes p and Kg mesons combined with random tracks and/or 
showers, these events will not interfere with each other or with resonances in the signal as 
they are not from a true D°. Additionally, Kgn° events populate a narrow region on the 
Dalitz plot in both signal and background. The corresponding amplitudes do not interfere 
with the other amplitudes that contribute to D° —> n + n~n° due to the long lifetime of the 
Kg. Therefore, the normalization of the Kg contribution component floats in the fit, but 
has no further role. The background shape is parameterized by a two dimensional cubic 
polynomial in (m 2 +7r _, m 2 +7r0 ) with terms representing p and Kg mesons. 


Figure [21 shows the Dalitz-plot distribution for the D° —> 7r + 7r _ 7r° candidates. Only con¬ 
tributions from p T 7r ± , p°7r° and K°n° are readily apparent. Modeling the background shape 
and correcting for efficiency as described above, we then parameterize the 7r + 7r _ 7r° Dalitz-plot 
distribution following the Breit-Wigner formalism using the unbinned likelihood method as 
described in Ref. jl5|,[l6j]. In a separate fit, we also parameterize the nn S-wave with the K- 
matrix formalism as described in Ref. 
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for the analysis of D + —* 7r + 7r - 7r + . We allow the 
normalization of the background contribution to float unconstrained in our fits. We consider 
seventeen resonant components, an 0 , p°n°, p Zf n ± , uin°, f 0 (980)n°, f 2 ( 1270)7t°, / 0 (1370)7r 0 , 
p(1450)°vr 0 , p(1450) =F 7r ± , / 0 (1500)7r°, p 3 (1690) 0 7r°, p(1700) 0 vr°, p( 1700) =f 7t ± , f 0 (1710)n°, as 
well as a non-resonant contribution. All interfere coherently and we fit for a complex co¬ 
efficient (amplitude and relative phase) for each resonance as well as for the non-resonant 
contribution. We describe the resonances with the standard parameters 0 . Lacking theo¬ 
retical guidance, the non-resonant contribution is modeled as a uniform distribution across 
the allowed phase space. 

This study is sensitive only to relative phases and amplitudes. The largest mode, p + n~, 
is assigned a zero phase and an amplitude of one. Since the choice of normalization, phase 
convention, and amplitude formalism may not always be identical for different experiments, 
fit fractions are reported in addition to amplitudes. The fit fraction is defined as the integral 
of a single component (resonant or non-resonant) over the Dalitz plot divided by the integral 
of the coherent sum of all components over the Dalitz plot 0 . The sum of the fit fractions 
for all components will in general not be unity because of interference. 

We use the full covariance matrix from the fits to determine the statistical errors on fit 
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FIG. 1: Distribution of a) M within 650 keV of our measured value of Q and b) Q within 44 MeV/c 2 
of our measured value of M for the process D° —> 7r + 7r _ 7r°. The candidates pass all selection criteria 
discussed in the text. The curves show the results of the fits. The vertical lines denote the signal 
region. The light (dark) shaded region indicates the signal (background) contribution. 


fractions to properly include the correlated components of the uncertainty on the amplitudes 
and phases. After each fit, the covariance matrix and final parameter values are used to 
generate sample parameter sets. The distributions of fit fractions from these parameter sets 
are then used to determine the Gaussian width and 95% confidence level (C.L.) upper limits. 

The results of our fits are presented in Table|T| Fit A includes the three p(770) resonances 
and an interfering non-resonant component. The non-resonant contribution is small and we 
do not find strong evidence for any other contributing resonances. 

We fit the D° and D° samples separately in Fit B1 and Fit B2, respectively. The violation 
of CP could manifest as distinct amplitudes and phases for D° and D° Dalitz plots. Since 
Fits B1 and B2 are consistent, there is no indication of CP violation. 

There has been significant interest in the properties of the mi S- wave due to the possibility 
of a low mass a meson and glueball of mass ~ 1.5 GeV/c 2 , each of which are beyond the 
qq' quark model [lij. A recent paper Hi highlights D decay as a preferred way to focus on 
these states because of the preponderance of A-waves in the initial and final state. Charged 
D mesons are observed to decay preferentially to the m r S- wave in the decay D + — > 7r + 7r - 7r + 


studied by FOCUS [11( and E791 [20]. We considered several possible contributions of a tttt 


S -wave component shown in Fig. 0 In Fit Cl we replace the flat non-resonant amplitude 
with a <j( 500) 0 contribution parameterized as a Breit-Wigner resonance. The contribution 
of a is very small in our fit. 

In Fit C2, we replace the <r(500) with the / 0 (980) where we use the mass and width deter¬ 
mined from the single Breit-Wigner parameterization of E791 [2l|. There are predictions fJ 
for D + —> / 0 (980)7t + (BF=2.8%) and D a —> / 0 (980)7r° (BF=0.06%). The prediction for the 
latter decay is extremely small because of the non-gg nature of / 0 (980) and strong final state 
interactions. We find a fit fraction of (0.010 ± 0.008)% for D° —> / 0 (980)7t°. Although sta¬ 
tistical errors limit the precision of our result, it is clear that this decay is highly suppressed 
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FIG. 2: a) The Dalitz-plot distribution for D° —> tt +/ k~tt 0 candidates. b)-d) Projections of the 
results of the fit described in the text to the 7r + 7r _ 7r° Dalitz plot showing both Fit D (line) and 
the data (points). The results of Fit A, Fit Cl, and Fit C2 are indistinguishable from Fit D. See 
text for details of the fits. 


in agreement with the model. 

The 7T7r S -wave has contributions from a number of overlapping resonances and there 
are several models which parameterize this wave from other reactions 0, 00- The 
complicated structure is much more amenable to a coupled channel formulation than Breit- 
Wigner models. In Fit D we consider a mr S -wave contribution following the A-matrix 
formalism of Au, Morgan, and Pennington 0 in addition to the resonant components of 
Fit A. As in Fit Cl and Fit C2, we exclude the non-resonant amplitude. Figure |2>d shows 
the three projections of the Fit D. These fits are very similar to those of FOCUS |i 11] for 
D + 


n Ti Ti 


~ -7T + 


The fit fraction of the 7T7 t S'-wave component in Fit Cl is much less than 
that of E791 for the analogous fit to D + —> 7r + 7r _ 7r + . For Fit D, our fit fraction, (0.9±0.7)%, 
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FIG. 3: Projections of the 7T7 t S'-wave parameterization for a) Fit A: flat non-resonant b) Fit Cl: 
Breit-Wigner <r(500) c) Fit C2: Breit-Wigner /q( 980) and d) Fit D: K-matrix. 


was very small when compared to the FOCUS value for D + decay, (56.00 ± 3.24 ± 2.08)%. 
The observed ratio of m t S -wave fit fractions in D + relative to D 0 is 60or > 36 @95% 
C.L. - somewhat larger than the tree-level estimate of (3\/2) 2 = 18. The factors in the 
relative amplitudes of 3 and are due to color suppression and isospin, respectively. 

The lack of evidence for a tttt S -wave in this analysis is interesting. The measured rates for 
the comparable reactions D + —* 7r + 7r°, (0.26±0.07)%, and D° —> 7r°7r°, (0.084±0.022)%, are 
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both given qualitatively by model calculations [8( which provide values of 0.19% and 0.11%. 
A quark model with final state interactions is used to fit parameters to a few charm decays 
and then predict many others. The same model predicts a ratio of branching fractions (BF) 
for D + —> / 0 (980)7t + and D° —> / 0 (980)7t° of 46.7. Comparing the results from Ref. j2lj | 
and our Fit C2 we measure this ratio to be 620^210 or > 340 @ 95% C.L.. Although the 
quantities measured in the experiments are not the same as what is calculated in the model, 
there is a strong likelihood that D three-body decay data can provide interesting tests of 
the microscopic nature of the scalar states. 

We calculate an integrated CP asymmetry across the Dalitz plot as described in Ref. Si] 


as the difference between the integral of the coherent sum of all amplitudes across the Dalitz 
plot for D° and D°, respectively, divided by the sum of the integral of the coherent sum 
of all amplitudes across the Dalitz plot for D° and D°, respectively, divided by the area 
of the Dalitz plot. We obtain Acp — O.Olj^io? ± 0.05, where the errors are statistical and 
systematic, respectively. 

We consider systematic uncertainties from experimental sources and from the decay model 
separately. Contributions to the experimental systematic uncertainties arise from our model 
of the background, the efficiency, the signal fraction and the event selection. Our general 
procedure is to change some aspect of Fit A or Fit D and interpret the change in the values 
of the amplitudes, phases, and fit fractions as an estimate of the systematic uncertainty. 
In Fit A and Fit D, we fix the coefficients of the background determined from a sideband 
region. To estimate the systematic uncertainty on this background shape we perform a fit 
with the coefficients allowed to float constrained by the covariance matrix of the background 
fit. We use the covariance matrix of the efficiency fit to estimate the systematic uncertainty 
due to the efficiency parameterization. We generate samples of efficiency parameters based 
on the covariance matrix and re-run the Dalitz-plot fit for each sample. We change selection 
criteria in the analysis to test whether our simulation properly models the efficiency. We vary 
the minimum 7T° daughter energy criteria, the cuts on Q and M, and the D* + minimum 
momentum fraction. We take the square root of the sample variance of the amplitudes, 
phases and fit fractions from the nominal result compared to the results in this series of fits 
as a measure of the experimental systematic uncertainty. 

We consider the uncertainty arising from the choice of the tttt S -wave model included in 
the fit. We interpret the variation in the p amplitudes, phases, and fit fractions in Fit A, 
Cl, C2 and D as a modeling systematic uncertainty. We add the experimental and model 
systematic uncertainty in quadrature to obtain the total systematic uncertainty reported in 
Table □ 

In conclusion, we have analyzed the resonant substructure of the decay D° —> 7r + 7r _ 7r° 
using the Dalitz-plot analysis technique. We observe the three p(770)7r resonant and a small 
non-resonant contribution. We find no evidence for a tttt S -wave contribution with either the 
Breit-Wigner or Ji-matrix parameterization. We find Acp which is the asymmetry between 
the D Q and D° distributions integrated over the entire Dalitz plot to be 0.011 q!o7 ± 0.05. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 



dation and the U.S. Department of Energy. 
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TABLE I: Fit A includes the three p( 770) resonances and a non-resonant (NR) contribution. We 
fit the D° and D° samples separately in Fit B1 and Fit B2, respectively. Fit Cl and Fit C2 allow 
a ct( 500) and /o(980) contribution, respectively, parameterized as a Breit-Wigner resonance [r| . 
The tt-k S-wscve contribution in Fit D is parameterized following the li-matrix formalism |TTi |. 



Amplitude 

Phase(°) 
Fit A 

Fit fraction(%) 

p + 

1. (fixed) 

0. (fixed) 

76.5±1.8±2.5 

p° 

0.56±0.02±0.03 

10±3±2 

23.9±1.8±2.1 

p~ 

0.65±0.03±0.02 

176±3±2 

32.3±2.1± 1.3 

NR 

1.03±0.17±0.12 

77±8±5 

2.7±0.9±0.2 




<6.4 @ 95% C.L. 



Fit B1 


P + 

1. (fixed) 

0. (fixed) 

76.6±2.5 

P° 

0.57T0.03 

104=4 

24.9±2.4 

P~ 

0.64T0.03 

1764=4 

31.0±2.8 

NR 

1.03T0.24 

724=11 

2.8T1.4 



Fit B2 


P + 

1. (fixed) 

0. (fixed) 

76.0±2.7 

P° 

0.55T0.04 

94=4 

22.5±2.7 

P~ 

0.67T0.04 

1774=4 

34.0T3.02 

NR 

1.03T0.24 

844=11 

2.7T1.4 



Fit Cl 


P + 

1. (fixed) 

0. (fixed) 

78.0±2.1 

P° 

0.56T0.02 

94=3 

24.4±1.9 

P~ 

0.66T0.03 

1764=3 

33.9±2.3 

cr(500) 

0.22T0.06 

3554=24 

0.08±0.08 




<0.21 @ 95% C.L. 



Fit C2 


P + 

1. (fixed) 

0. (fixed) 

78.3±1.8 

P° 

0.56T0.02 

104=3 

24.9±1.9 

P~ 

0.66T0.03 

1784=3 

33.4±2.1 

/o(980) 

0.074T0.025 

325 ±23 

0.010±0.008 




<0.026 @ 95% C.L. 



Fit D 


P + 

1. (fixed) 

0. (fixed) 

76.3±1.9±2.5 

P° 

0.57±0.03±0.03 

10±3±2 

24.4±2.0±2.1 

P~ 

0.67±0.03±0.02 

178±3±2.0 

34.5±2.4± 1.3 

A-matrix 0.70±0.20±0.12 

2±14±5 

0.9±0.7±0.2 




<1.9 @ 95% C.L. 
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